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The  objective  of  this  paper  is  to  study  the  spreading  of  a  submerged  two- 
phase  jet.  A  plunging  liquid  jet  impacting  on  a  still  liquid  pool  may 
carryunder  bubbles  by  entraining  the  surrounding  gas.  For  low  liquid 
jet  turbulence,  the  measured  average  bubble  size  was  about  200  m 
[Bonetto  and  Lahey,  1992].  For  low  void  fraction,  one  may  assume  that 
the  liquid  velocity  field  is  not  disturbed  by  the  presence  of  the  bubbles. 

This  allows  the  use  of  known  solutions  for  the  velocity  field  in  a  planar 
liquid  jet.  That  is  the  liquid  flow  may  be  solved  using  Eulerian 
coordinates  and  neglecting  the  effect  of  the  bubbles.  The  momentum 
equation  for  the  bubbles  has  been  written  in  Lagrangian  coordinates. 

The  resulting  equations  are  solved  for  bubble  trajectories  in  the  known 
liquid  velocity  field. 

INTRODUCTION 

The  gas  entrapment  by  a  plunging  liquid  jet  and  the  resultant  two-phase 
jet  dispersion  occur  in  many  problems  of  practical  interest.  The 
ecological  balance  in  lakes  and  oceans  is  dependent  on  the  amount  of 
dissolved  oxygen.  In  addition,  the  air  carryunder  that  occurs  at  most 
hydraulic  structures  in  rivers  is  primarily  responsible  for  the  large 
air/water  mass  transfer  that  is  associated  with  these  structures  [Avery 
and  Novak,  1978].  Also,  the  absorption  of  CO2  by  the  oceans  appears  to 

play  an  important  role  (as  a  sink  term)  in  reducing  the  magnitude  of  the 
"greenhouse  effect".  The  absorption  of  greenhouse  gases  has  been 
hypothesized  to  be  highly  dependent  upon  the  air  carryunder  that  occurs 
due  to  breaking  waves. 
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Plunging  jets  are  also  used  for  dissolving  gas  into  liquid.  Some 
chemical  reactors  use  a  plunging  jet  to  entrain  bubbles,  the  reaction  rate 
is  increased  because  of  the  increase  in  interfacial  area  and  the 
enhancement  in  the  mass  transfer  characteristics  produced  by  the 
turbulence  in  the  two-phase  jet. 

The  problem  under  study  is  a  planar  liquid  jet  impacting  on  a 
liquid  pool  surrounded  by  a  still  gas.  Lezzi  and  Prosperetti  [1991]  and 
Bonetto  et  al  [1992]  have  recently  proposed  mechanisms  for  the  gas 
entrainement.  In  this  work  the  dispersion  of  a  two-phase  jet  is  studied. 
A  dilute  dispersion  of  bubbles  in  the  liquid  is  assumed.  The  liquid 
velocity  field  is  assumed  to  not  be  affected  by  the  presence  of  the 
bubbles.  The  momentum  balance  equations  are  presented  in  the 
Momentum  Balance  Section.  Then,  the  Tolmien's  solution  for  the  mean 
velocity  of  a  turbulent  liquid  jet  is  derived.  Finally,  the  integrated 
bubble  trajectories  and  other  derived  quantities  of  interest,  such  as  the 
maximum  penetration  depth  of  the  bubble,  the  total  time  that  the  bubbles 
are  under  water  and  their  emergence  position  at  the  pool's  surface  are 
presented. 

MOMENTUM  BALANCE 


Newton's  law  for  the  gas  phase  takes  the  form 


P  V 
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where  g,  V^,  Ug  are  the  particle  density,  volume  and  velocity, 

D  Ug 

respectively.  — ^ — is  the  material  derivative  of  the  particle  velocity 

Dt  Ao«a«ni(Ui  for 

following  the  particle  and  F  is  the  applied  force. 

!  IVM  tAS 

It  is  convenient  to  divide  the  force,  F,  into  three  parts,  one  being  the ! 

inviscid  flow  force  (Re=oo)  and  the  others  corresponding  to  viscous  flows  | — ^ - 

Fy  and  the  buoyancy  force  Fb-  i 

F  =  El  +  Ev  +  Eb  (2)  L '  "  *  ♦ 
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i)  Inviscid  flow  force:  The  inviscid  flow  force  can  be  further  separated 
into  two  parts,  one  corresponding  to  an  irrotational  flow  and  the  other 
corresponding  to  a  rotational  flow. 


„(irrot)  (rot) 

El  =  El  +  El 


(3) 


For  spherical  bubbles  in  inviscid  irrotational  flows  the  force  El  assuming 
that  the  liquid  phase  pressure  gradient  is  impressed  on  the  bubbles,  is 
given  by: 


Dt 


(4) 


D.u^ 

where  U/  is  the  liquid  velocity,  — —  is  the  material  derivative  following 

Dt 

the  liquid  and  Cym  is  the  virtual  mass  coefficient  (for  spheres  Cvin=l/2). 
The  term  is  zero  because  in  this  study  the  liquid  flow  was  steady, 
thus  the  liquid  acceleration  is  given  by. 


=  (u.  •  V)u^ 


A  general  expression  has  been  derived  for  the  lateral  force  on  a  sphere  in 

IVuJR. 

weakly  rotational  flows  (Vu^  ^  « 1).  Auton  et  al  [1988]  and 

Drew  &  Lahey  [1987,1990]  show  that, 

=  Cj^(Ug -u^)x  Vxu£  (5) 


where  Cl  is  the  lift  coefficient  which,  for  spheres  is,  Cl=1/2. 

ii)  Viscous  Flows  (moderate  Re):  The  effect  of  the  viscous  stress  will  now 
be  considered.  A  particle  moving  in  a  steady  flow  has  a  non-steady 
component  of  the  drag  force  due  to  vortex  shedding.  The  viscous  force  is 
then: 
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Fv(t)  =  FD+F^(t) 


(6) 


E0  =  ip,CpA^IUrlUr  (7) 

Where  Ed  is  the  steady  drag  force  component  parallel  to  the  relative 
velocity  and  FjL(t)  is  the  non  steady  force  component  perpendicular  to  the 

relative  motion.  For  moderate  Re  numbers  (Re  =  100)  the  drag  force  is 
caused  by  the  pressure  field  asymmetry  in  the  flow  around  the  particle. 
The  vortex  shedding  produces  an  oscillating  pressure  distribution  around 
the  bubble  and  a  non-steady  force  component.  Graham  [1980]  showed 
that  if  the  Strouhal  number  is  S=Rb/(lui-l  T)  <  0.3,  Fj_(t)  is  less  than  10  % 
of  the  value  of  Ed,  where  T  is  the  period  of  the  vortex  shedding.  In  our 
simulation  the  previous  condition  is  satisfied  and  hence  Ei(t)  is  neglected. 
That  is, 

Ev  =  Ed  (8) 

iii)  Buoyancy  force:  The  buoyancy  force  is  due  to  the  difference  in 
density  between  the  bubbles  and  the  liquid, 

EB  =  (P^-pg)Vbg  (9) 


Finally,  we  note  that. 


Dt 


(10) 


where  x  is  the  position  of  the  particle  at  time  t.  The  momentum  balance 
equation  is  obtained  by  combining  Eqs.  (1)  to  (5),  (8)  and  (9).  It  is 
convenient  to  write  Eq.  (10)  and  the  momentum  balance  in  component 
form  for  a  two-dimensional  steady-state  flow: 
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5  P/bCo  I  u;  -  Ug  I  ( u<  -  Ug )  +  ( -  Pg ) Vfa  g 
+  Cj^pgVbUjX(Vxu<)  (11c) 

TOLMIEN  SOLUTION 

Let  US  now  consider  the  solution  for  a  spreading  turbulent  single-phase 
liquid  jet.  The  jet  has  a  uniform  velocity  profile  where  it  impacts  the 
liquid  pool  at  z=0.  After  the  impact,  a  boundary  layer  develops  at  the 
edge  of  the  jet.  Figure- 1  shows  the  different  jet  regions.  Region- 1  is  the 
inviscid  core.  Inside  this  region  the  velocity  is  the  inlet  velocity,  V. 
Region-2  is  where  the  jet  is  approximately  self-similar.  At  larger 
distances  from  the  impacting  point  the  jet  is  closer  to  a  true  self-similar 
condition.  Tolmien's  calculation  will  be  followed  to  obtain  [Hinze,  1975] 
the  mean  velocity  values  in  the  self-similar  region.  It  is  assumed  that  the 
liquid  emanates  from  an  infinitely  narrow  cross  section.  Also  it  is 
assumed  that  the  pressure  inside  the  jet  is  uniform  and  equal  to  the 
surrounding  pressure.  Then  the  momentum  crossing  a  plane 
perpendicular  to  the  direction  of  the  flow  is  constant: 

pu^dy  =  constant  (12) 

•^0 

Assuming  the  velocity  profile  is  self-similar. 
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(13) 


U  Y 

=  then 

u  z 


rb/z' 


^0 


idy 

(-^r  “7  =constant 
z' 


where  z'  is  the  axial  distance  measured  from  the  infinitely  narrow  cross 
section.  Combining  eqs.  12  and  13  we  obtain  the  relationship  between 
the  center-line  axial  velocity  uz(y=0)  and  z*  is  as: 

where  c  is  a  constant.  The  velocity  at  an  arbitrary  point  is  given  by: 
u^=^f(n)  (15) 


It  is  useful  to  define  the  stream  function  as: 


dw  d\if 

u  =  — ^  u  =  — — 
z  9y  y  9z 


(16) 


then 

\|f  =  |u  dy  =  cVz'F('n)  , where, 

Z 

F(Ti)  =  Jf(Ti)d'n 


(17) 


The  momentum  balance  in  the  z  direction  for  the  control  volume  shown 
in  Figure-1  is: 


y 

r 


U  U  +T— 

z  y  9z' 


u^dy  +  d^z'^ 


'^0u  ^ 
_ z 


8y 


=  0 


(18) 


where  d  is  the  mixing  length  constant  and,  according  to  mixing  length 
theory,  the  turbulent  shear  stress  was  assumed  to  be. 
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(19) 


T  =  -pd 


rau  f 

_ z 

ay , 


Using  Eqs.  (16)  and  (17)  we  obtain: 

2d2(F"(Tl))2  =  F(Tl)F'(tl)  (20) 

Introducing  the  additional  change  of  variables, 

(p  =  —  , where,  a  =:^(2d^)  (21) 

a 

yields, 

(F"(<p))2  =  F(<p)F'(<p)  (22) 

The  appropriate  boundary  conditions  at  the  centerplane  (p=0  are  that  the 
lateral  velocity  Uy  equals  zero  and  the  axial  velocity  Uz  is  Uz(y=0),  thus, 

F((p  =  0)  =  0  ,  F'(cp  =  0)  =  l  (23) 

Equation  (22)  has  been  numerically  integrated  with  the  boundary 
conditions  Eq.  (23).  In  order  to  compute  the  velocities  the  value  of  the 
parameter  has  to  be  known.  The  value  obtained  experimentally  by 
Forthmarm  [Abramovich,  1963]  was  used,  that  is,  a=0.11. 

Note  that  this  calculation  corresponds  to  a  jet  that  is  self-similar 
everywhere  and  that  it  emerges  from  an  infinitely  narrow  slit.  Actually, 
the  jet  emerges  from  a  finite  slit  with  half  width  equal  to  h  and  constant 
velocity  V.  In  order  to  reconcile  these  differences  an  effective  origin  is 
computed  for  the  Tolmien  jet,  and  it  is  required  that  its  momentum  agree 
with  the  actual  jet  momentum  at  impact.  The  momentum  flux 
conservation  condition  is: 
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(24) 


P/ 


f  u^dy  =  p^v2h 
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The  distance  from  the  infinitely  narrow  slit,  z\  and  the  distance  from  the 
finite  slit,  z,  are  related  by, 

z'=z+s  ,where,  s=0.41h/a 

A  final  modification  has  to  be  made  to  the  self-similar  solution  to  account 
for  the  inviscid  core  region.  For  the  infinitely  narrow  solution  the 
velocities  goes  to  infinity  we  get  closer  to  the  slit.  In  the  finite  slit,  the 
velocity  can  not  exceed  the  inlet  velocity  V,  thus  wherever  the  velocity 
exceeds  V,  it  is  reset  to  V.  Figs.  2  and  3  show  the  lateral  and  the  axial 
velocity  Uy  and  xiz  as  a  function  of  the  lateral  coordinate  y  for  different 

axial  positions  z.  Note  that  at  the  inlet  the  velocity  profile  is  uniform  and 
the  jet  is  spreading  as  we  increase  the  axial  position  z.  Figs.  4  to  8  show 
the  various  spatial  derivatives. 


DRAG  LAW 


To  this  point,  the  analysis  is  valid  for  both  bubbles  and  solid  particles. 
However,  there  is  a  difference  between  the  drag  that  the  liquid  exerts  on  a 
bubble  compared  to  that  for  a  solid  particle. 


The  so-called  "dirty  water"  Wallis  drag  correlation  as  been  used  to 
compute  the  drag  on  bubbles: 


Cd- 


6.3 

ReO.385 


(25) 


In  contrast  the  Clift,  Grace  and  Weber  correlation  is  generally  used  to 
compete  the  drag  in  liquid/solid  flows 


^(l  +  0.15ReP«*n 
Re'  ' 


0.42 


H-4.2510'^Re'‘'-^^] 


(26) 
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where  Re=2  lurl/\)^ 

RESULTS 

Before  presenting  the  results  is  convenient  to  make  the  following  change 
of  variables: 

=  y  cos(0)  +  z  sin(0) 

=  y  sin(0)  -  z  cos(0) 

Figure  9  shows  the  trajectories  of  tliree  bubbles  with  initial  positions  (a) 
(y=-h/2,z=Rb),  (b)  (y=0,z=Rb)  and  (c)  (y=h/2,z=R5).  Bubble  (a)  is 
initially  at  the  center  plane,  while  the  other  two  are  initially  at  one  side  or 
the  other  of  the  jet  centerplane  halfway  to  the  edge  of  the  jet.  For  the 
case  shown  in  Figure  9,  the  liquid  planar  jet  is  impacting  the  pool  at  an 
angle  of  0=45  degrees  with  respect  to  the  vertical.  The  uniform  velocity 
of  the  jet  at  the  impacting  point  is  V=5m/s,  the  jet  width  is  2h=4.16  mm 
and  the  radius  of  the  bubble  is  Rb=2mm.  For  the  base  case  we  used  the 

coefficients  of  the  lift  force  and  the  virtual  mass  force  computed 
analytically  assuming  that  the  liquid  was  inviscid,  that  is,  Cvm=l/2  and 
Cl=1/2. 

A  qualitative  description  of  the  flow  follows.  We  consider  first  the 
bubble  that  is  initially  at  the  center  of  the  jet.  In  the  inviscid  core  region, 
all  derivatives  of  the  liquid  velocity  field  are  zero,  then  the  only  forces 
acting  acting  on  the  bubbles  are  drag  and  buoyancy.  The  drag  force  acts 
on  the  bubble  in  the  same  direction  that  the  liquid  flow.  In  the  inviscid 
core  the  lateral  velocity,  Uy^  is  zero  and  the  axial  velocity,  u^,  is  equal  to 

V.  Thus,  the  drag  force  in  the  inviscid  core  tends  to  move  the  bubble 
parallel  to  the  jet  and  at  a  gas  velocity  near  V.  The  buoyancy  force  acts  in 
the  upward  vertical  direction.  That  is,  buoyancy  tends  to  cause  the  bubble 
to  move  towards  the  free  surface.  In  the  self-similar  region,  the  lift  force 
keeps  the  bubble  near  the  centerplane,  the  virtual  mass  reduces  the 
acceleration  that  the  bubble  experiences  for  a  given  force. 

We  take  the  previous  solution  as  the  base  solution  in  parameter  space. 
Then  the  effect  of  the  different  parameters  are  explored  one  at  a  time  and 
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compared  with  the  base  solution. 

MAXIMUM  DEPTH 

One  of  the  important  parameters  is  the  maximum  depth  that  a  given 
bubble  reaches.  Figure  10  shows  the  maximum  depth  as  a  function  of  the 
inlet  liquid  velocity  V  for  three  different  inlet  positions.  The  increase  of 
the  maximum  depth  with  the  velocity  is  easily  seen.  This  is  expected 
because  when  we  increase  the  velocity  by  a  factor  AV,  the  lift  force  that 
keeps  the  bubble  near  the  center  plane  increases  a  factor  (AV)^  while  the 
buoyancy  force  that  tends  to  move  the  bubble  toward  the  free  surface 
remains  the  same.  Notice  that  the  depth  reached  by  bubble  with  initial 
conditions  in  the  lower  half  of  the  jet  have  similar  trajectories  compared 
to  the  bubbles  initially  at  the  center  plane.  This  is  again  due  to  the  fact 
that  the  bubbles  that  remain  near  the  centerplane  longer  are  the  ones  that 
reach  greater  depths.  A  typical  trajectory  of  a  bubble  that  is  at  the  lower 
half  of  the  jet  initially  is  as  follows:  buoyancy  tends  to  move  it  up  (i.e. 
toward  the  center  plane),  however  as  the  bubble  gets  closer  to  the  center 
plane,  the  higher  liquid  velocities  there  produce  a  larger  drag  and  the 
bubble  flows  toward  greater  depths. 

EMERGENCE  DISTANCE 

Another  quantity  of  importance  is  the  distance  from  the  jet  impact 
position  to  the  place  where  the  bubble  reemerges  at  the  free  surface. 
Figure  1 1  shows  the  distance  of  emergence,  s,  as  a  function  of  the  jet  inlet 
velocity  V  for  three  different  initial  bubble  positions.  Again  the  5 
increases  with  the  inlet  velocity  V  because  of  the  effect  of  the  lift  force  . 
It  is  clearly  seen  in  Figure  11  that  the  bubbles  that  are  initially  in  the 
lower  half  of  the  jet  (yi<0)  have  a  longer  emergence  distance  s  than  the 

bubbles  initially  in  the  upper  half  of  the  jet  (y>0).  Therefore  s  increases 
as  a  function  of  the  initial  bubble  position  yo  (-h<yo<h).  Thus,  the  larger 
the  initial  position  yo  (-h<yo<h),  the  larger  the  emergence  distance  s, 
except  for  yo  very  near  -h.  For  yo  values  near  -h,  an  interesting 
phenomenon  occurs.  Figure  9  shows  two  typical  trajectories  with  yoj>yo2 

both  in  the  neighborhood  of  yo  =-h.  The  buoyancy  force  acts  over  the 
bubble  with  initial  position  yoj  producing  a  lateral  velocity  U/y  >0  which 
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is  greater  than  the  lateral  velocity  of  bubble  number  2,  Uy^  >Uy2  >0. 

Then,  bubble- 1  crosses  the  jet  centerplane  sooner  than  the  bubble -2. 
Thus,  bubble  2  has  a  larger  emergence  distance  s. 

Notice  that  the  two  trajectories  cross  each  other  .  An  Eulerian  approach 
would  give  some  kind  of  average  behavior  at  this  point.  This  result  does 
not  contradict  the  fact  that  a  set  of  Ordinary  Differential  Equations 
(ODE)  has  an  unique  solution.  The  trajectories  are  non-crossing  in  full 
phase  space.  The  dimension  of  the  space  is  equal  to  the  number  of 
independent  equations  (i.e.,  four).  Thus  the  coordinates  in  phase  space 
are  y,z,Uy,Uz.  The  trajectories  presented  in  the  y-z  plane  can  be  seen  as 

projections  of  the  trajectories  in  full  phase  space  on  the  y-z  plane.  Note 
that  the  xi-x2  plane  is  identical  to  the  y-z  plane.  This  effect  can  not  occur 

using  partial  differential  equations  resulting  from  the  Eulerian  approach. 
The  corresponding  solution  would  effectively  average  the  different 
velocities  at  trajectory  crossings,  and  would  result  in  a  single  velocity  at 
each  point. 

SENSITIVITY  OF  THE  SOLUTION  TO  THE  PROBLEM 
PARAMETERS 


We  wish  to  study  the  effect  of  the  various  parameters  on  the  resulting 
bubble  trajectories.  What  follows  is  a  parametric  study  of  the  solution  of 
eq  (1 1).  In  the  following  figures  the  dashed  lines  correspond  to  the  base 
solution  trajectories  from  the  previous  section.  The  solid  lines  show  the 
effect  of  changing  the  corresponding  parameters.  Again  curves  (a),  (b) 
and  (c)  correspond  to  the  initial  conditions  described  in  Fig.  9. 

In  Figs,  13  and  14,  the  lift  coefficient  was  varied,  with  Cl=0.25  (for 
a  base  case  solution  of  Cl=0.5).  For  the  trajectories  (a)  and  (b)  the 
reduction  in  the  lift  coefficient  produces  a  reduction  in  the  ability  of  the 
bubble  to  remain  near  the  center  plane.  The  y  component  of  the  lift  force 
is: 


F("o‘)  =  r  (u.  -u  ) 
ly  L"  «  g/ 


3u, 


_ y 


dy  dz 


<0,  if,  y>0 
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Figure  8  shows  the  derivatives  as  a  function  of  y.  This  factor  is  always 
negative  for  y>0.  Also  the  axial  component  of  the  relative  velocity  is 
positive.  Thus  for  y>0  the  lift  force  tries  to  push  the  bubblea  to  the 
centerline.  At  the  center  plane  the  axial  liquid  velocity  has  a  maximum 
and  therefore  the  drag  acts  to  move  the  bubble  to  greater  depths. 
Obviously,  there  is  a  strong  dependence  of  the  emergence  point,  s,  on  the 
value  of  the  lift  coefficient. 


There  is  also  a  lift  force  component  in  the  axial  direction: 


F(rot)  ^ 
I 

z 


V 


3a^  'I 

_ J_ 

9z 
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For  bubbles  in  the  upper  half  of  the  jet,  the  lateral  bubble  velocity 
component  is  in  general  greater  than  the  lateral  liquid  velocity 
component.  Thus,  the  lift  force  component  in  the  axial  direction  is 
negative.  For  trajectory  (c)  the  negative  lift  force  component  in  the  z 
direction  is  reduced  and  the  emergence  position,  s,  is  increased.  In  Fig  14 
the  lift  force  has  been  completely  eliminated  (Cl=0)  and  we  find  that  the 

emergence  distance  of  the  trajectories,  s,  are  fairly  close  together.  Hence 
we  find  that  a  larger  lift  coefficient  results  in  a  spreading  of  the 
emergence  distance  for  the  various  trajectories. 

Figs.  14  and  15  correspond  to  varying  the  virtual  mass  coefficient 
from  the  base  case  value  of  =  0.5,  to  Cym  =0-4  and  Cym  =1.0.  We 

note  that  a  larger  virtual  mass  coefficient  results  in  less  acceleration  of  the 
bubble,  so  that  an  increase  in  the  virtual  mass  coefficient  produces  an 
increase  in  the  ability  of  the  bubble  to  follow  the  liquid  flow. 

Figs.  16  and  17  corresponds  to  the  effect  of  the  jet  inclination  angle, 
0.  Fig.  16  shows  a  jet  impacting  at  an  angle  of  40  degrees  (base  case 
=45°).  The  maximum  depth  dm  is  seen  to  be  sensitive  to  the  angle  .  On 

the  contrary,  the  location  of  the  emergent  bubble  is  seen  to  be  somewhat 
insensitive  to  the  angle  .  Fig.  18  shows  a  jet  impacting  at  an  angle  =60 
degrees.  For  this  case  both  the  maximum  depth,  dm.  and  the  distance  of 

emergence,  s,  are  seen  to  be  strongly  affected  by  the  impacting  angle. 
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SUMMARY  AND  CONCLUSIONS 


A  dilute  bubbly  jet  has  been  analyzed  using  a  Eulerian-Lagrangian 
approach.  It  is  assumed  that  the  presence  of  the  bubbles  does  not  affect 
the  liquid  flow.  Under  this  assumption  the  liquid  flow  is  decoupled  from 
the  bubble  motion,  and  is  given  by  the  Tollmien  self-similar  jet  solution 
modified  to  account  for  the  inviscid  core.  That  is,  it  was  not  necessary  to 
solve  for  the  liquid  flow  field.  The  momentum  balance  for  the  gas  phase 
has  been  rewritten  as  a  system  of  four  ODEs.  The  independent  variable  is 
time  and  the  four  dependent  variables  are  the  bubble  position  (y,  z)  and 
the  bubble  velocity  components  (uy,  u^).  The  system  of  ODEs  was 

integrated  using  a  fifth  order  Runge  Kutta  scheme.  The  initial  conditions 
were:  y(t=0)=yo,  z(t=0)=0,Uy(t=0)=0,Uy(t=0)=V. 

Based  on  the  computed  trajectories,  three  quantities  of  practical 
interest  were  estimated:  The  maximum  depth  d^j,  the  bubble  emergence 

distance,  s,  and  the  total  time  that  the  bubble  is  immersed  in  the  pool,  t. 
One  may  be  interested  in  knowing  the  aeration  produced  by  a  breaking 
wave  in  the  sea.  The  maximum  depth  provides  information  about  the 
amount  of  the  aerated  layer.  The  time  under  water  gives  information 
about  the  amount  of  gas  that  is  dissolved.  However,  in  order  to  compute 
the  amount  of  dissolved  gas,  an  additional  model  for  the  diffusion  process 
needs  to  be  included.  Finally,  the  distance  of  emergence  gives 

information  about  the  horizontal  extent  of  the  aeration  produced  by  the 
breaking  of  a  single  wave.  The  results  show  that  d^,  t  and  s  increase  with 

the  jet  velocity,  with  faster  increase  for  low  values  of  V.  Also  the  bubbles 
initially  at  the  lower  half  of  the  jet  have  larger  dm,  t  and  s  than  the 
bubbles  initially  at  the  upper  half  of  the  jet  for  the  same  velocity. 

The  relative  velocity  in  the  axial  direction,  Ug^-u^^  is  negative  inside 

the  jet.  Thus  the  lateral  component  of  the  lift  force  tends  to  move  the 
bubble  toward  the  jet  centerplane.  The  axial  component  of  the  lift  force 
is  I  '*<?ative  if  the  relative  velocity  is  positive  and  this  force  tends  to  retard 
the  bubble  motion. 

The  effect  of  virtual  mass  force  is  not  as  straightforward  to  analyze. 
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The  reason  is  that  the  terms 


and  ^  change  sign  inside  the  liquid 
ay  dz 


jet.  This  can  be  easily  seen  in  Figs.  4  and  7. 


We  conclude  that  both  the  virtual  mass  and  the  lift  forces  play  an 
important  role  in  the  spreading  of  a  two-phase  jet. 
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1 :  inviscid  core 
2:  self-similar  region 


Figure  1 :  The  inviscid  and  self  similar  regions  of  a  planar  jet 


Figure  2:  Lateral  velocity  u.  as  a 
different  axial  positions  z 
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z=0.2m 

(c)  z=0.3m 

(d)  z=0.4m 
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function  of  the  lateral  position  y  for 


(a)  z=0.1in 

(b)  z=0.2m 

(c)  z=0.3m 

(d)  z=0.4m 


Figure  3:  Axial  liquid  velocity  as  a  fimction  of  the  lateral  position  y  for 
different  axial  positions  z 
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Figure  5:  — —  as  a  function  of  the  lateral  position  y  for  different  axial 
dz 


positions  z 
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Figure  7:  as  a  function  of  the  lateral  position  y  for  different  axial 

oz 

positions  z 


(a)  z=0.1m 

(b)  z=0.2m 

(c)  z=0.3m 

(d)  z=0.4m 


Figure  8:  Multiplier  for  the  lift  coefficient  as  a  function  of  the  lateral  position 
y  for  different  axial  positions  z 


Rb=0.002m 


Figure  9:  Bubble  trajectories  for  the  base  solution  and  different  initial 
positions. 

(a)  y(t=0)=-hy2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 
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Figiire  10:  Maximum  depth  dm  as  a  function  of  the  jet  velocity  for 
different  initial  positions: 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 


s(m) 


Figure  11:  Emergence  distance  s  as  a  function  of  the  jet  velocity  for 
different  initial  positions: 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=K))=+h/2 


Figure  12:  Emergence  time  t  as  a  function  of  the  jet  velocity  for  different 
initial  positions: 

(a)  y(t=0)=-h/2 
y(t=0)=0 
(c)  y(t=0)=+h/2 
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0=45° 

V=5ni/s 

Rb=0.002m 


Figure  13:  Bubble  trajectories  for  Cl=0.25  and  different  initial  positions 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 

The  dashed  lines  are  the  base  solutions 


Rb=0.002m 


Figure  15:  Bubble  trajectories  for  Cvin=0-4  and  different  initial  positions 
(a)  y(t=0)=-h/2 
y(t=0)=0 
(c)  y(t=0)=+h/2 

The  dashed  lines  are  the  base  solutions  (0^^  =  0.5) 


Rb=0.002m 


Figure  16:  Bubble  trajectories  for  Cvm=l-0  different  initial  positions 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 

The  dashed  lines  are  the  base  solutions  (C^  =  0.5) 


Rb=0.002in 


Figure  17:  Bubble  trajectories  for  0=40°  and  different  initial  positions 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 

The  dashed  lines  are  the  base  solutions  (6=45°) 


Rt)=0.(X)2m 


Figure  18:  Bubble  trajectories  for  0=60°  and  different  initial  positions 

(a)  y(t=0)=-h/2 

(b)  y(t=0)=0 

(c)  y(t=0)=+h/2 

The  dashed  lines  are  the  base  solutions  (0=45°) 


